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^ The random intensity of noise approach to ID Laval-DubruUe-Nazarenko model is used to describe 

Lagrangian acceleration of a fluid particle in developed turbulence. Intensities of noises entering 
nonlinear Langevin equation are assumed to depend on random velocity fluctuations in an exponen- 

1 tial way. The conditional acceleration PDF, variance, and mean are found to be in a good qualitative 
, agreement with the recent high-Re Lagrangian experimental data. 
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Q Phcnomcnological approaches 0, were used to describe Lagrangian acceleration of a fluid particle in developed 
D turbulent flow within the framework of Langevin type equation; see also Recent one-dimensional (ID) stochastic 
^ , particle models and some refinements @ were reviewed in • Some toy models of developed turbulence suffer from 
^ I ^ • the lack of physical interpretation deduced from turbulence dynamics Isl . 

[ Recently 0, 0] we have shown that the ID Laval-DubruUe-Nazarenko (LDN) type toy model [lolflll | of the accelera- 
• tion evolution with the model turbulent viscosity and coupled delta-correlated Gaussian multiplicative and additive 
I ; ' noises is in a good agreement with the recent high-precision Lagrangian experimental data on acceleration statistics; 
d Rx ~ 690, the normalized acceleration range is [—60,60], Kolmogorov scale is resolved ll^ ll^ ll ^ . The longstanding 
^ ' Heisenberg-Yaglom scaling, (a^) = aou^^^i'~^^'^L~^/'^ , was confirmed experimentally |12| to a very high accuracy, 
I , for about seven orders of magnitude in the acceleration variance, or two orders of the rms velocity u, at R\ > 500. 
' ^ Long-time correlations and the occurrence of very large fluctuations at small scales dominate the motion of a fluid 
particle, and this leads to a new dynamical picture of turbulence ll Sj, Il6l. We focus on modeling the acceleration 
statistics conditional on velocity fluctuations u presented recently in 

The original 3D and ID LDN models were formulated both in the Lagrangian and Eulerian frameworks for small- 
scale velocity increments in time and space respectively. They are based on a stochastic kind of Batchelor-Proudman 
rapid distortion theory approach to the 3D Navier-Stokes equation |lO|. and thus have a deductive support from 
, turbulence dynamics. The random intensity of noise (RIN) approach 0, 0] provides an extension of the above ID 
' model in the limit of small time scale r for which Lagrangian velocity increments are proportional to r: u(t + t) — 
I u{t) = Ta{t). The main idea of RIN approach is simply to account for the recently established two well separated 
T-H ■ Lagrangian autocorrelation time scales for the velocity increments [15| and assume that certain model parameters, 
[ such as intensities of the noises, fluctuate at the long time scale due to Lagrangian velocity u. A simple ID model 
. can shed some light to properties of 3D LDN model of Lagrangian dynamics. 

We give only a brief sketch of the ID LDN model and refer the reader to l^i llo| for more details; see also [l^ . 
This toy model can also be viewed as a passive scalar in a compressible ID flow [l^- We use probability density 
function (PDF) obtained as a stationary solution of the Fokker-Planck equation associated to the Langevin equation 
for the component of Lagrangian acceleration a{t): da/dt = {£, — vtk^)a + a±_ ^^l- ■'-'^ 3D LDN model, ^(i) is related 
' to the velocity derivative tensor and cr±(i) describes a forcing of small scales by large scales via the energy cascade 
' mechanism. In ID LDN model, these are approximated by external Gaussian white noises, 

O 

^ ■ m) = 0, mm) = ^OSit ~ f), {a^{t)) = O, {a^{t)a^{t')) = 2aS{t - t'), = 2A<5(< - t'), (1) 



that is partially justified by DNS The acceleration PDF can be calculated exactly 0, with the result 



H ' P{a) = Cexp[-i^tk^/D + F{c) + F{~c)]{Da^ -2\a+af''^{2Bka + i^tk^)-^^^''/^ , (2) 



for constant parameters. Here, we have denoted ft — \/ i^q + B'^a^/k^, C is normalization constant, 



^(c) = TT-TvT 7 + - ^")« + c{Dvlk^ -f C2i^t))], c = -i^Da^X^, (3) 

IC2IJ c CiC2{c — JJa + A) 
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ci = B^{4\^+Ac\'^-3DaX - cDa) + D'^ic+X)iy^k'^ , = ^J B'^{2X^ + 2cA - Da)B + ZJV^fc^. (4) 

Without loss of generality one can put, in a numerical study, k — \ and the additive noise intensity q = 1 by 
rescaling of the multiplicative noise intensity > 0, the turbulent viscosity parameter B > 0, the kinematic viscosity 

> 0, and the cross correlation A 0, and make a fit of P{a) to the experimental data. The particular cases B — 
and i/Q = at A = were studied in detail in . Nonzero A is responsible for an asymmetry of the PDF |(5J and in 3D 
picture corresponds to a correlation between stretching and vorticity (the energ y ca scade) ; in the Eulerian framework, 
{{6u)^) is proportional to A, in accord to a kind of Karman-Howarth equation 

Since the experimental unconditional and conditional distributions, Pexp(a) and Pexp(a|w) at w = 0. were found to 
be approximately of the same stretched exponential form revealing strong Lagrangian intermittency |l4| we use the 
result of our fit Q of the PDF to Poxp(a) measured with 3% relative uncertainty for \a\ < 10. This implies 
the following rescaled parameter set: D = 1.100, B = 0.155, vq = 2.910, A = -0.005 (fc = 1, a = 1, C = 3.230). 
A fit to the conditional distribution Pcxp(a|0) would yield a different set of values of the parameters. The used fit 
is however justified on a qualitative level. We assume that the parameters a and A entering ^ depend on the 
amplitude of (normalized) Lagrangian velocity fluctuations u, while D, B, and are taken to be fixed at the fitted 
values (fc = 1). An exponential form of a{u) has been proposed in J] and was found to be relevant from both the 
(K62) phenomenological and experimental points of view. Particularly, such a form leads to log-normal RIN model 
when u is independent Gaussian distributed with zero mean (PDF g{u)), and yields the acceleration PDF whose 
low-probability tails are in agreement with experiments 0, 0; the marg inal PDF is P,n(a) = J^^ P{a\u)g{u)du. 
Remarkably, this form is found to provide appreciable increase of the conditional acceleration variance (a^jw) with 
increasing |u| 7] that meets the experimental data |l4j . Guided by these observations the simplest choice is to try an 
exponential dependence for A(u). Particularly, in the present description we use X{u) = -0.005e3l"l and a{u) = el"l. 
This allows us to deal with the acceleration statistics which exhibits a strongly non-Gaussian character in both the 
unconditional and conditional cases. 

The conditional PDF is given by Q treated in the form P{a\a{u), A(u)). The constant C in (|2J) is calculated for 
each value of u. For the normalized velocity fluctuations values u = 0, 0.25, 0.50, 0.75, 1.00, 1.19 the distributions 
P{a\a{u), X{u)) and mean accelerations (a|u)/(a^|0)^/^ are shown in Fig-s.^andl^J One observes a good qualitative 
correspondence of the obtained conditional PDFs (Fig. ^ shifted for clarity) with the experimental curves (Fig. 6a 
in 101). Both the variance and skewness of P{a\u) increase for bigger While the increase of the variance (related 
to the increase of a{u)) is readily seen, the increase of the skewness (related to the increase of X{u)) results in a 
rather small change of the shape of distribution despite the fact that A varies by about two order of magnitude (from 
|A|=0.005 to 0.18). Such a change can be however readily seen in a plot of the contribution to fourth order moment, 
a'*P(a) (see Fig. 4 in |^). The obtained conditional mean acceleration (a|u)/(a^|0)^/^ plotted in Fig. |21is also in a 
good qualitative agreement with the experimental dependence (a|u)/(a^)^/^ (Fig. 6b in [III). The mean acceleration 
is zero for a symmetrical distribution and for homogeneous isotropic turbulence. One observes a rather small relative 
increase of the mean acceleration for bigger |w| that eventually reflects a coupling of the acceleration to large scales 
of the studied flow [TolIT^. We note that the experimental (a|u)/(a^)^/^ exhibits small asymmetry with respect to 
u — > —u. 
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In summary, the presented ID model is shown to capture main features of the observed conditional acceleration 
statistics. 
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